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Abstract. A nonhomogeneous system of linear recurrence equations can be recognized by an automaton J?l over a one-letter 
■ alphabet A = {z|. Conversely, the automaton generates precisely this nonhomogeneous system of linear recurrence equations. 

We present the solutions of these systems and apply the z-transform to these solutions to obtain their series representation. Finally, 
we show some results that simplify the series representation of the z-transform of these solutions. We consider single systems as 
well as the composition of two systems. 
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. 1. Introduction 

Consider a finite automaton J[ over a one-letter alphabet A - {z}, with weights in a semiring R. Since A = {z}, an 
arbitrary word w e A* has the form w = z ■ Z ■ • -Z = Z n , for some n e N. By definition, the weight of a word w - z" equals 
the sum of the weights of all paths with label w — z" 12 16| . But since every transition reads the letter z, the weight of the 
word z" equals the sum of the weights of all paths of length n. Note that, in this kind of automaton, we are not interested 
in accepting/rejecting words over some alphabet, but rather in counting paths of length n, and keeping track of the weights 
of such paths. The idea of using automata as counting devices has been used recently with applications in combinatorics 
and difference equations BHUEI, but goes back to classical automata theory Q. 
jy-j \ Now suppose we are interested in computing the weight of all paths of length n (equivalently, all paths with label z") in 

an automaton 3K, starting at a specific state p. Then we would look at all paths starting at p, compute the weight of every 
path of length n, and add up the results. Denote this function by f p , and notice that for every n e N, f p (ri) gives the sum of 
the weights of all paths of length n starting at state p. Hence, in this kind of automaton, every state p generates a function 

Let 3\ be a weighted £-state automaton over A = {z} and an arbitrary semiring R, and let /i,/2, . . . ,fk 6 R N be the 
functions generated by states p\,pi, . . . ,pk, respectively. Then it can be shown [1 , 4| that these functions satisfy the 
following system of homogeneous linear recurrence equations with constant coefficients: 



X 
S3 



/i(n + 1) = fln/i(n) + anfzin) + ... + aufk(n) 
f 2 (n + 1) = a 2 \f\(n) + anfiin) + ... + a 2 kfk(n) 

fk(n + 1) = a k \f\{n) + a k 2fi{n) + ...+ a kk fk(n) 

where a,y 6 R is the weight of the transition p,- — > pj. Conversely, given this system of linear recurrence equations, the 
weighted finite automaton recognizing it is precisely ZR.. 

This result can be generalized to consider the case when the system is nonhomogeneous and has variable coefficients 
IT]. That is, there is an automaton J\ over A = {z} and R that generates the following system of linear recurrence equations: 

fi(n + 1) = a n {n)f\{n) + a l2 (n)f 2 (n) + ...+ a lk (n)f k (n) + gi(n) 
f 2 (n + 1) = a 2 \(n)f\(ri) + a 22 (n)f 2 (n) + ...+ a 2k (n)fk(n) + g 2 (n) 



f k {n + 1) = a k i(n)fi(n) + a k2 (n)f 2 (n) + ... + a kk (n)f k (n) + g k (n) 

l 
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where fi,gi,atj € R M . And conversely, given this system of linear recurrence equations, the automaton recognizing it is 
precisely 3K. The structure of this automaton is different from the structure of the automaton recognizing the homoge- 
neous system. The difference accounts for the fact that the nonhomogeneous system contains two types of functions; the 
functions gi are the input whereas the functions /) are the output. A state in J\ stores the current value of the function 
it represents. If we identify a state with the function it generates, then ay e R H is the variable weight of the transition 
fi < — fj and, for every i, there is an additional transition gi — > /;, with weight 1 e R. This kind of automaton generalizes 
the one recognizing the homogeneous system as it recognizes the nonhomogeneous system, as well as the homogeneous 
one. 

a n (n) 
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(Mn)) 
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, g(n) = 




, and A(n) = 


Jk(n), 




,gk(n), 





a 2 \(n) 



flki(n) a k2 (n) 



diking 
aikin) 



a kk {n). 



Then we can rewrite the above system 



as f(n + 1) = A(n)f{n)+g(n), where f,ge Rf xl and A e Rf xk . 



Sometimes we consider the (simpler) case when the system 

and 



has constant coefficients. In that case, we are dealing with the system /(« + 1) = Af(n) + g(n), where /, g e Rf xl 

A 6 Rkxk- 

We will first consider the solutions of these systems, and then move on to consider the solutions of composing two 
systems. We will then apply the z-transform to these solutions to obtain their series representation. Finally, we will show 
a result that presents a simplified series representation of their z-transform. 



2. Solutions of Nonhomogeneous Systems of Linear Recurrence Equations 

From now on, assume for simplicity that the initial conditions of the systems above are 0. That is, we assume that 
/(0) = 0. Then it can be shown [1 1 that the solution of the system with constant coefficients is given by 



(2.1) f(n + 1) = A" * g(n), 

where * denotes discrete convolution: for functions a(n) and B(n), {a * B)(n) = Y!l=o a ( n ~ ■ /?(/)■ I n order to give the 
solution of the system with variable coefficients, we need a definition. 

Definition 2.1. The symbol A {k) (n) is defined by A w (n) = A{n)A(n - 1) • ■ • A(n - (k - 1)), for k = 1,2, ...,«+ 1. This 
definition can be extended to allow k — by letting A (0 \n) = I, the identity matrix. (Note that k represents the number of 
matrices being multiplied.) 

Using this definition, it can be shown [ 1] that the solution of the system with variable coefficients is 



(2.2) /(« + l) = [A w (?) *£(«)],=„. 

In the expression [A (n) (f) * g(n)] t= „, the exponent n is the variable of the convolution, used to expand the expression 
A (n) (f) * g(ri), while the argument t stays fixed until after the convolution has been expanded, and then is substituted by 
t = n. That is, [A«(f) * g{n)] l=n = [A«(t) ■ g(0)]t=n + [A^(t) ■ s(l)],=„ + . . . + [A«(t) • g(n - 1)] (= „ + [A<°>(f) • g(n)] t=n = 
A(n) • ■ -A(l)g(0) + Mn) ■ ■ ■ A(2)g(l) + . . . + A(n)g(n - 1) + g(n). 

Now suppose we have two k by k systems /(n + 1) = Af(n)+g(n) and g(n + 1) = Bg(n)+h(n), with constant coefficients. 
We can compose these two first-order systems to obtain a second-order system. The solution of this system gives / in 
terms of h. It can be shown [ lj that its solution is 



(2.3) f(n + 2) = A" *[B" * h(n)]. 

Similarly, we could consider two k by k systems with variable coefficients: /(« + 1) = A(n)f(ri) + g(n) and g(n + 1) = 
B(n)g(n) + h(n). The result of composing these two systems yields / in terms of h. It can be shown [ 1| that its solution is 



(2.4) f(n + 2) = [A (n) (f) * [B (n \s) * h(n)] s=n ] t=n+1 . 

Example 2.1. Let f(n + 1) = (n + l)/(n) + g(n) and g(n + 1) = ng(n) + h(n) be systems over R = N, where /(0) = 0, 
g(0) = and h(n) — n + 1. Note that A(n) — n + 1 and B(n) = n. Suppose we want to find /(4). First, direct computation 
shows 



And by using this, we obtain 



g(l)=l,g(2) = 3,g(3) = 9. 
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/(l) = 0,/(2)=l,/(3) = 6,/(4) = 33. 
We can check this result by using formula (12.4b . which allows us to compute / directly from h: 

f(4) = f(2 + 2) = [A«(f) * [B ( "\s) * h(n)] s=n ] l=n+l \ n=2 

= [A< 2 >(f)[B <0) (s) * h(0)] s=0 + A^(t)[B^(s) * h(l)] s=l + A<°)(f)[B (2) (^) * h(2)] s=2 ] t=3 

= [A (2) (f)[B (0) (0)/j(0)] + A (1) (f)[B (I) (l)/;(0) + B (0) (l)/i(l)] + A (0) (f)[B (2) (2)/!(0) + B (I) (2)/;(1) + B (0) (2)/i(2)]], =3 
= A (2) (3)[B (0) (0)/!(0)] + A (1) (3)[B (1) (l)/i(0) + B (0) (l)/i(l)] + A (0) (3)[B (2) (2)/i(0) + B (1) (2)/;(l) + B {Q \2)h{2)} 

= 12[l]+4[l+2] + l[2+4 + 3] =33 

Formulas d2.lt and (12. 2t will allow us to compute the z-transform of / when / is given by a system of linear recurrence 
equations. Formulas (12. 3t and ( 12. 4t will allow us to compute the z-transform of / when / is given by the composition of 
two systems of linear recurrence equations. 

3. The z-Transform of the Solution of a System 

The z-transform is a useful tool when working with functions defined as sequences, as in our case. We now define the 
z-transform. 

Definition 3.1. Given a function 

/ = (/«)), /(l ),..., /(«),... )eR N , 
we define the z-transform of / to be the series 



/(0)/z° + /(D/z 1 + . . . + f{n)lz n + ••• = £ f(n)/z". 
We denote the z-transform of / by z(f). 



n=0 



Suppose now that / is the k by 1 vector of functions / 



h 
\fk> 



, defined by / = 



f(/i(0),/i(l ),..., fi («),... )\ 
(/ 2 (0),/ 2 (l ),..., fiin),...) 



e R 



kxV 



We can apply the z-transform to every component, thus representing / as 



Alternatively, we can write 



TT-l DO 

as 2j„=o 



fi(n)\ 
fiin) 

Mn), 



l(A(0),A(l),...,A(n),...)J 
/i(0)/z° + /i(D/z I + ... + /i(»)/z" + . ^ 

/2(0)/z°+/ 2 (l)/z I + ...+/2(«)/z" + . 

MO)/z° + Ml)/z 1 + ... + Mn)/z n + . 
/z", via an isomorphism. Note that this last 



fzr=o/i(")/^ 

Z: =0 f2(n)/z n 

expression is equal to Y*7=o fi n )/z n - Hence we can say that, in the general k dimensional case, the image of / e /?rL under 
the z-transform is Yj7=o fi n )/z"- 

In the constant case we are given a k by k system 

f(n + l)=Af(n)+g(n); 

in the variable case we are given a A: by A: system 

f(n + l)=A(n)f(n)+g(ri). 

The solutions to these systems are, respectively, 

fin + 1) = A" * g(ri), 

and 

f{n + \) = [A (n \t)*g{n)] t=n . 
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These two expressions we have for fin + 1) (one for the constant and one for the variable case) give the coefficients in the 
z-transform of /. 

Since /(0) = 0, we obtain in both cases that zif) = £~ /(")/z" = /(«) + ZZo f(» + D/z" +1 = 2T=o /(« + D/z" +1 ■ In 
the constant case, 

OO 

(3.1) z(/) = X [A "*S (n)]/z " +1 ' 
and in the variable case, 

oo 

(3.2) zif) = £[A W « * S(n)Wz" +1 . 

n=Q 

Using the definition of convolution, it is possible to rewrite the above series in a simpler fashion. In the constant case, 

z(f) = Z: =Q lA n *g(n)]/z n+1 

= [A * g(0)]/z l + [A 1 * g(l)]/z 2 + [A 2 * g(2)]/z 3 + [A 3 * g(3)]/z 4 + . . . 

= [AV0)]/z' 
+[A 1 g(0)+A°g(l)]/z 2 
+[A 2 #)+AV(D+AV2)]/z 3 
+[A 3 g(0) + A 2 g(l) + A l g{2) + A°g(3)]/z 4 + ... 

= [A°/z° + A 1 /? 1 + A 2 lz 2 + . . -]g(0)/z l 
+[A°/z° + A'/z 1 + A 2 /z 2 + . . .]g(l)/z 2 
+[A°/z° + A'/z 1 + A 2 /z 2 + . . .]g(2)/z 3 + . . . 



= [Z2 =Q A m /z m ]g(0)/z l 



Z% =0 A m /z" 

oo 

m=0 y 



2: =0 A m /z m 



g(D/z 2 
#(2)/z 3 + 



Z: = o[Z% = oA m /z m }g(n)/z"+ l 



Notice that the series Yjm=oA m /z m is part of the coefficient of l/z" + ', but it does not depend on n. We can rename this 
series as S a = Yjm=oA m /z m , and then write the series for zif), in terms of g, as zif) = Yj™=o S Agin)/z n+1 - This proves the 
following theorem. 

Theorem 3.1. Let fin + 1) = A/(n) + gin) be a k by k system of nonhomogeneous linear recurrence equations with 
constant coefficients, with initial conditions equal to 0. Then the z-transform of f is 

OO 

(3.3) zif) = J]s A gin)/z n+l , 

where S A = Y2^A m lz m . 
In the variable case, 

zif) = ZZo\-A (n) it)*gin)] t =n/z n+1 

= [A<°>« * gi0)] t=0 /z l + [A^it) * g(l)], = i/z 2 + [A^it) * g(2)] f=2 /z 3 + [A< 3 >(f) * #(3)] (=3 /z 4 + . . . 

= [A(°)(f)g(0)] (=0 /z 1 
+[A< 1 >(%(0) + A<°>(f)s(l)L=i/z 2 
+[A (2 \t)giO) + A^iOgil) + A<°)(f)g(2)], =2 /z 3 
+[A< 3 >(f)g(0) + A^it)gi\) + A^iOgil) + A^it)gi3)] 1= i/z 4 + ... 

= [A(°)(0)g(0)]/z 1 
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+[A( 1 >(l)g(0) + A(°>(l)g(l)]/z 2 

+[A (2) (2)g(0) + A^i2)gil) + A ( -°\2)g(2)]/z i 

+[A( 3 >(3)g(0) + A< 2 >(3)g(l) + AW(3)g(2) + A (Q H3)g(3)]/z 4 + ... 

= [A<W(0)/z° + A< 1 >(l)/z 1 + A< 2 >(2)/z 2 + . . .]g(0)/z l 
+[A<°>(l)/z° + A < - l \2)/z l + A®(3)/z 2 + . . .]g(l)/z 2 
+[A<°>(2)/z° + A^QVz 1 + A®(4)/z 2 + . . .]g(2)/z 3 + . . . 



2: = oA (m) (m+l)/z" 
2: = oA (ffl) (m + 2)/z" 



S(D/z 2 
g(2)/z 3 + ... 



= JZa [JZ* A<m> ( m + n)lz m \g{n)lz^ 

Note that the series Zm=o^ (m) ( m + n)/z m , which is part of the coefficient of 1 /z" +1 , is a function of n. By letting Sa(i) = 
2™ =0 A (m) (m + n)/z m , we can write the series for z(/), in terms of g, as zif) = Yj™ = q S A(n)gin) / z" +1 ■ This shows the 
following theorem. 

Theorem 3.2. Let fin + 1) = A(n)/(n) + g(n) be a k by k system of nonhomogeneous linear recurrence equations with 
variable coefficients, with initial conditions equal to 0. Then the z-transform of f is 

oo 

(3.4) z(/) = ^Sa(»)s(»)/z" +1 , 

where S A (n) = \Z2 =() A (m \m + n)/f\ 

4. The z- Transform of the Solution of the Composition of two Systems 
Suppose / is defined by 

f(n + 1) = A f{n) + gin), where g(n + 1) = Bg(n) + /<«), 

in the constant case; 

fin + 1) = A(n)/(«) + g(n), where g(n + 1) = Bin)gin) + hin), 
in the variable case. In both cases, we can solve for / in terms of h. The solutions are 

fin + 2) = A n * [B" * hin)], 

in the constant case, and 

fin + 2) = [A ( ">(0 * [B<»>( S ) * hin)] s=n ] t=n+u 
in the variable case. Again, notice that these expressions give the coefficients in the z-transform of /. 

We assume that /(0) = and g(0) = 0. Thus, /(l) = 0. We obtain in both cases that zif) = X~=o/(»)/z B = 
/(0) + /(l) + Zr =0 /(» + 2)/z" +2 = 2r = o /(» + 2)/z" +2 . In the constant case, 

OO 

(4.1) zif) = Y} A " * * M«)]]/z" +2 , 

and in the variable case, 



(4.2) zif) = £[A W (0 * [B'»V) * ^(n)]. != „] ( = n+ i/z" 



n ^ * [u yt) f "V»7Jj=«J(=«+l/'" 

n=0 

Just as we did before, it is possible to write these series as simpler expressions. In the constant case, 



zif) = XZo[A n *[B n *hin)]]/z n+2 

= [A * (B° * /z(0))]/z 2 + [A 1 * (B 1 * /z(l))]/z 3 + [A 2 * (B 2 * /z(2))]/z 4 + . . . 

= [A°iB°hiO))]/z 2 
+[A 1 iB°hiO)) + A°iB l hiO)) + A°(B°/z(l))]/z 3 

+[A 2 (fi°/z(0)) + A 1 (5^(0)) + A\B°hil)) + A°(fi 2 /z(0)) + A Q (B l h(l)) + A°(fi°/z(2))] /z 4 + . . . 
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: A°fi°/z(0)/z 2 

+A 1 B°h(0)/z 3 + A Q B 1 h(0)/z 3 + A°B°h(l)/z 3 

+A 2 B°h(0)/z 4 + A 1 B 1 h(0)/z 4 + A 1 B°h(l)/z 4 + A°B 2 h(0)z 4 + A°B l h(l)/z 4 + A°B°h(2)/z 4 + . . 

A°B°h(0)/z 2 + A 1 B Q h(0)/z 3 + A°B l h(0)/z 3 + A 2 B°h(0)/z 4 + A l B x h{Q)lz 4 + A°B 2 h(0)/z 4 + . . 
+A°B°h(l)/z 3 +A 1 B°h(l)/z 4 + A°B l h(l)/z 4 +A 2 B°h(l)/z 5 +A 1 B l h(l)/z 5 + A°B 2 h(l)/z 5 + . 
+A°B°h(2)/z 4 + A l B°h(2)/z 5 + A°B 1 h(2)/z 5 + A 2 B°h(2)/z 6 +A l B 1 h(2)/z 6 + A°B 2 h(2)/z 6 + . 

■■ [A°B°]h(0)/z 2 + [A l B°/z l +A°B 1 /z l ]h(0)/z 2 + [A 2 B° / z 2 + A 1 B l / z 2 + A B 2 / z 2 ]h(0) / z 2 + ... 
+[A°B°]h(l)/z 3 + [A 1 B°/z l +A°B 1 /z l ]h(l)/z 3 + [A 2 B° / z 2 + A 1 B 1 / z 2 + A B 2 / z 2 ]h(l) / z 3 + . . 
+[A°B°]h(2)/z 4 + [A l B Q lz x +A°B 1 /z l ]h(2)/z 4 + [A 2 B Q / z 2 + A 1 B 1 / z 2 + A B 2 / z 2 ]h(2) / z 4 + . . 

■■ [(A * B Q )/z°]h(0)/z 2 + [(A 1 * B l )/z l ]h(0)/z 2 + [(A 2 * B 2 ) / z 2 ]h(Q) / z 2 + ... 
+[(A° * B°)/z°]h(l)/z 3 + [(A 1 * B l )lz x ]h{\)lz 3 + [(A 2 * B 2 ) / z 2 ]h(l) / z 3 + ... 
+[(A° * B°)/z°]h(2)/z 4 + [(A 1 * B 1 )/z 1 ]h(2)/z 4 + [(A 2 * B 2 ) / z 2 ]h(2) / z 4 + ... 



= ['L2 = o(A m *B m )/f t ]h(0)/z 2 



TZ=o(A m * B m )lf 
YZ =Q (A m * B m )lf 



KD/z 3 
K2)/z 4 ■ 



ji+2 



= JZa [Z2= Q (A m * B m )/f] h(n)/z n 

The series H~ =0 (A m * B m )/z m is part of the coefficient of 1 /z" +2 and it does not depend on n. Hence, we can let S a,b = 
Zm=o( A ' n * B m )/z m and write the series for z(/), in terms of h, as z(f) = S A , B h{n) I f +2 . We have shown the following 
theorem. 

Theorem 4.1. Let fin + 1) = A f(ri) + g(ri) and gin + 1) = B fin) + h(ri) be two kbyk systems of nonhomogeneous linear 
recurrence equations with constant coefficients, with initial conditions equal to 0. Then the z-transform of f is 

oo 

(4.3) zif) = J]s A , B hin)/z" +2 , 



n=0 



where S A , B = ^(A" 1 * B m )/f\ 
In the variable case, 

zif) = Zr = o[A w « * [B^is) * hin)] s=n ] t=n+l /z n+2 

= [A«(0 * [B^is) * hin)] s=n ] t=n+l \ n=Q /z 2 
+[A { "\t) * [B^\s) * Kn)] s=n ] t=n+1 \ n=l /z 3 
+[A M (t) * [Bf*\s) * h(n)] s=n ] t=n+l \ n=2 /z 4 + . 



= [A ( - \t)[B^is)*h{0)] s=0 ] t=1 /z 



2 

+[A^it)[B^\s) * h(0)] s=0 + A<°>(f)[B (1) (s) * /z(l)]. s=1 ], =2 /z 3 

+[A^ 2 \t)[B^is) * h(0)] s=0 + A^\f)[BP-\s) * h(l)l s=1 + A<P\t)[B^ 2 \s) * K2)] s=2 ] t=3 /z 4 + ... 

= [A<°>(f)[5 (0) (0)/i(0)]], =1 /z 2 
+[A (1) (/)[fi (0) (0)/z(0)] + A (0) (f)[B (1) (l)/i(0)] + A (0) (f)[B (0) (l)/z(l)]],=2/z 3 
+[A (2) (f)[fi (0) (0)/z(0)] + A (1) (f)[B (1) (l)/z(0)] + A (1) (f)[B (0) (l)/z(l)] 
+A (0) (f)[£ (2) (2)/z(0)] + A (0) (f)[fi (1) (2)Ml)] + A (0) (f)[B (0) (2)/z(2)]], =3 /z 4 + . . . 

= [A(°)(l)[BW(0)/z(0)]]/z 2 
+[A (1) (2)[B (0) (0)/z(0)] + A (0) (2)[B (1) (l)/z(0)] + A (0) (2)[B (0) (l)/z(l)]]/z 3 
+[A (2) (3)[B (0) (0)/i(0)] + A (1) (3)[B (1) (l)/z(0)] + A (1) (3)[B (0) (l)/z(l)] 
+A (0) (3)[B (2) (2)/z(0)] + A (0) (3)[B (1) (2)/z(l)] + A (0) (3)[B (0) (2)/i(2)]]/z 4 + . . . 
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= A (O \l)B (O \0)h(0)/z 2 
+A (1) (2)fi (0) (0)/z(0)/z 3 + A (0) (2)fi (1) (l)/z(0)/z 3 + A (0) (2)fi (0) (l)/z(l)/z 3 
+A (2) (3)fi (0) (0)/z(0)/z 4 + A (1) (3)5 (1) (l)/z(0)/z 4 + A (1) (3)5 (0) (l)/i(l)/z 
+A (0 >(3)B (2 \2)h(0)/z 4 + A(°'(3)B (1) (2)/z(l)/z 4 + A< »(3)8 (0) (2)/i(2)/z 4 + . . . 

= A (0) (l)B (0) (0)/z(0)/z 2 + A (1) (2)B (0) (0)/z(0)/z 3 + A (0) (2)B (1) (l)/i(0)/z 3 
+A (2) (3)fi (0) (0)/z(0)/z 4 + A (1) (3)fi (1) (l)/z(0)/z 4 + A (0) (3)fi (2) (2)/z(0)/z 4 + . . . 
+A (0) (2)fl (0) (l)/z(l)/z 3 + A (1) (3)fi (0) (l)/z(l)/z 4 + A (0) (3)fi (1) (2)/z(l)/z 4 
+A (2) (4)5 (0) (l)/i(l)/z 5 + A (1) (4)5 (1) (2)/z(l)/z 5 + A (0) (4)5 (2) (3)/z 2 /z(l)/z 5 + . . . 
+A (0) (3)fi (0) (2)/z(2)/z 4 + A (1) (4)fi (0) (2)/z(2)/z 5 + A (0) (4)fi (1) (3)/i(2)/z 5 
+A( 2 '(5)B (0) (2)/z(2)/z 6 + A (r >(5)B^(3)h(2)/z 6 + A^ \5)B (2 \4)h(2)/z 6 + ... 

= [A <0 >(l)B^(0)/z°]h(0)/z 2 + [A <1) (2)B <0 \0)/z 1 +A(°»(2)B< 1 )(l)/z 1 ]/i(0)/z 2 
+[A (2) (3)fi (0) (0)/z 2 + A (1) (3)B (1) (l)/z 2 + A (0) (3)5 (2) (2)/z 2 ]/i(0)/z 2 + . . . 
+[A^(2)B (0 \l)/z°]h(l)/z 3 + [A^(3)B (0 \l)/z l + A<°>(3)fl (1) (2)/z 1 ]fc(l)/z 3 
+[A (2) (4)B (0) (l)/z 2 + A (1) (4)B (1) (2)/z 2 + A (0) (4)B (2) (3)/z 2 ]/i(l)/z 3 + . . . 
+[A^(3)B (0) (2)/z°]/i(2)/z 4 + [A^(4)B (0 \2)/z l +4(°'(4)B< 1 »(3)/z 1 ]/i(2)/z 4 
+[A (2) (5)B (0) (2)/z 2 + A (1) (5)B (1) (3)/z 2 + A (0) (5)B (2) (4)/z 2 ]/z(2)/z 4 + . . . 

= [AW(f) * B (m) (m)/z°] t=1 | m=0 /*(0)/z 2 + [A^(f) * B^\m)/z l ] t=2 \ m=1 h(0)/z 2 
+[A ( "»(() * B {m \m)lz 2 ] t=i \ m=2 h{Q)lz 2 + ... 

+[A^\t + 1) * B {m \m + l)/z°] f= i| m=0 Ml)/z 3 + [A w (( + 1) * B {m \m + I) I z l ] t=2 \ m= yh(l) I z* 
+[A«(f + 1) * B (m \m + l)/z 2 ] t= 3\ m=2 h(l)/z 3 + ... 

+[A<"% + 2) * B (m \m + 2)/z°] f=1 | m=0 /i(2)/z 4 + [A< m )(f + 2) * B«(m + 2)/z 1 ] f=2 | m=1 /i(2)/z 4 
+[A W (» + 2) * B< m >(m + 2)/z 2 ] t=3 | m=2 /z(2)/z 4 + . . . 



[Z^)[>1 W (0 * B (ffl) (m)],=m + i/z ffl ] K0)/z 2 
+ 
+ 



2: =0 [A (ffl) (f + 1) * B (m) (m + l)] (=m+1 /z" 
2 ; : =0 [A (m) (f + 2) * Z?( m )(m + 2)] f=m+I /z" 



KD/z 3 
h(2)/z 4 + . . 



2T = o [LZ= [A (m \t + n) * flW(m + »)] (=m+1 /z m ] h(n)li 



n+2 



The series 2m=o[ A(m) ( f + ") * B (m) (m + n)],= m+ i/z m , which is part of the coefficient of l/z" +2 , is a function of n. We 
can let S A , B (n) = 2m =0 [ A(m) ( f + «) * + nj] t=m+ i/f, and thus write the series for z(f), in terms of h, as z(/) = 

ZZ7=o S A,B(n)h(n) / z" +2 ■ We have proved the following theorem. 

Theorem 4.2. Let f(n + 1) = A(n)f(n) + g(n) and g(n + 1) = B(n)f(n) + h(n) be two kbyk systems of nonhomogeneous 
linear recurrence equations with variable coefficients, with initial conditions equal to 0. Then the z-transform of f is 



(4.4) z(f) = J]s A , B (n)h(n)/z n+2 , 

where S A , B (n) = ^[A^f + ») * B (m \m + n)], =m+l /z m . 
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